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We derive the formulae for the cross section of J/tp production in high energy pA and AA 
collisions taking into account the gluon saturation/color glass condensate effects. We then 
perform the numerical calculations of the corresponding nuclear modification factors and find 
a good agreement between our calculations and the experimental data on .1/^ production in 
pA collisions. We also observe that cold nuclear modification effects alone cannot describe the 
data on J/ip production in AA collisions. Additional final state suppression (at RHIC) and 
enhancement (at LHC) mechanisms are required to explain the experimental observations. 



I. INTRODUCTION 



The goal of this paper is to provide an improved analysis of the gluon saturation effects on the 
color singlet mechanism of J/ip production in dA and AA collisions at RHIC and LHC. In our 
recent publications [IHS] we argued that a mechanism responsible for J/^ production in central 
nuclear collisions is different from the one in pp collisions. This is because the symmetry properties 
of J/ip under the parity and charge conjugation transformations dictate that there must be an odd 
number of gluons attached to the bound c and c quarks. At the lowest order in strong couphng 
as there are three gluons attached. In pp collisions, two of those gluons have their external ends 
attached to the valence quarks of the colliding protons whereas the third one is emitted by the 
cc dipole. On the other hand in central pA collisions the parametrically enhanced contribution in 
the quasi-classical regime - which is controlled by a large parameter a'^A^^^ ~ 1 - originates 
from the diagrams where one of the gluons is attached to the proton's valence quark whereas the 
remaining two are attached to the valence quarks inside two different nucleons of the nucleus. 
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Obviously, such contribution breaks the perturbative QCD factorization aheady at the leading 
order in as- 

In [T] we assumed that the cc pair propagates through the nucleus in the color octet state and 
becomes color singlet only after the last interaction with the nucleus. In this paper we drop this 
assumption by taking into account a possibility that the cc pair converts from the color octet to 
the color singlet state already inside the nucleus. In the large Nc approximation further color 
conversions of the cc state are suppressed and thus can be neglected. Therefore, in this case the cc 
experiences the last inelastic interaction inside the nucleus after which it rescatters only elastically. 
As a result, the last inelastic interaction does not exponentiate with the rest of the scatterings 
and - as we will show - automatically selects an odd number of inelastic scatterings as required 
by the parity of J/ip- This is different from our approach in [IH3] where we had to select the odd 
number of inelastic scatterings in the scattering amplitude. Additionally, we give a more accurate 
treatment of J/ip wave function with parameters taken from a fit to the exclusive J/ip production 
in deep inelastic scattering. 

Our paper is structured as follows. In Sec. |ll]we derive the cross section for J/^ production 



in pA collisions; our main result is given by Eq. (15). In Sec. Ill we propose a generalization of 
this result to the A A collisions. The derived cross section is given by (17), (18) and satisfies the 
constraints imposed by the symmetry of the J/V' wave function. The results in Sec. |ll] and Sec. Ill 
are derived in the quasi-classical approximation, i.e. assuming that the coherence length for J/ijj 
production is much larger than the nuclear radius, but neglecting the low-x evolution. In Sec. 11^ 



we derive expression for the scattering amplitude (26) that includes the low-x evolution and thus 
gives a dependence on energy and rapidity. Sec. |V]is dedicated to the description of the numerical 
calculations performed with different models for the dipole scattering amplitudes. Our main results 



are exhibited in Figs. 3j4 We discuss them and conclude in Sec. VI 



II. PRODUCTION OF J/iP IN PA COLLISIONS 

The cross section for J/tp production in pA collisions can be written in the factorized form 

In order to set normalizations for '^"^"'^^^^'^^ it is convenient to compare the gA scattering process in 
([1]) with that of jA where there is a well developed phenomenology. Start with 7-proton scattering 
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where 



da. 



7P— s> J/f/jp 

~dt 



with 



A, 



7P— J/V'P 



(x,A) 



d^be 



-iA-b 



167r 



dz 



A. 



7p— ^ J/ipp\ 



d^r 
in 



rj/M'^i [l-5(x,r,b)] 



and t is given in terms of the momentum transfer hy t = — A^. Call 



(2) 



(3) 



(4) 



where 



2 

$^(r,z) = 776— {mlKo{mcr)(f)Tir, z) - \z'^ + (1 - 2)^] mcKi{mcr)dr4>T{r, z)] (5) 

O TT 



with [3 m 



iiTC?", -2) = — z) exp 



(6) 



and where Nt = 1.23, = 6.5 GeV-^ |8]. 

Except for a factor of z(l — z) in ([g]) our notation, and choice of J/il) wave function exactly 
matches that of Ref. [9]. Because ([T]) is a collinear factorized expression the gluon projectile on the 
right hand side of ([T]) is on-shell and so only transverse polarizations appear. We have taken the 
photon in ^ also on-shell so that the relationship between the photon and gluon induced processes 
will involve only a normalization change in ^ and a change of the 1 — S factor in ^ . 
We can get ([2]) in a more convenient form by using ([3]) and integrating over A. Thus 
da^ 



'-yA^J/tpA' 



l\z I ^,^^(r,z)l\z' I ^^;ir',z')[l-S*ir')][l-S{r)] (7) 



where we have suppressed the energy and impact parameter dependence in the 1 — S factors in 
The S factors are given, in the McLerran-Venugopalan model by 



S{r) = exp 



Cf^Qs^ 2 



exp 



(8) 



and the cross section in Q allows nuclear breakup but is elastic at the dipole-nucleon scattering 
level. Qs in ^ is the gluon saturation momentum with impact parameter dependence again 
suppressed. 

The main change necessary to convert ([7|) to a cross section for gA — )• J/iIjX is the way the 
cc dipole scatters off nucleons in the nucleus. In Q the scatterings are purely elastic, and such 
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FIG. 1: Smaple diagram contributing to the gA — > J/tp process. The point of the last inelastic interaction 
is signaled explicitly at the longitudinal coordinate ^. 

scatterings are dominant in the large-A'^c limit because the quantum numbers of the 7 and the 
J/ip are the same. In gA collisions the cc pair emerging from the gluon is in the adjoint color 
representation. The cc forming the J/V' is, of course, a color singlet. In the large- A'^c approximation 
there is a particular dipole-nucleon inelastic collision which converts the adjoint representation to 
a color singlet. This inelastic interaction is at the longitudinal coordinate ^, starting from the front 
of the nucleus, in Fig. [l} Later interactions, occurring after the cc pair is in a singlet state, are 
purely elastic in order to keep the singlet intact. Earlier interactions, occurring while the cc is 
in the adjoint representation, may be either elastic, occurring off a single c or c in the amplitude 
or complex conjugate amplitude, or inelastic involving the c in both the amplitude and complex 
conjugate amplitude or involving the c in both the amplitude and complex conjugate amplitude. 
Sample interactions are illustrated in Fig. [TJ 
The interaction at ^ gives the factor 

Qlr-r' 



4T{b) 

The interactions occurring before ^ give the factor 



-di. (9) 



while those occurring after give 



In going from 7^4 — ^ J/ipA' to gA — ^ J/ij^X the [1 — S*{r')] [1 — 5(r)] factor in ^ gets replaced 
by the product of the factors in ([9|-(11). 

In addition there is a color factor. In 7 induced J/'ip production there is a factor of Nc in the 
amplitude and a factor of N^, in the complex conjugate amplitude. This is the factor of Nc explicit 
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FIG. 2: Lowest order process in gluon induced J production. Color indices are indicated explicitly. 

in ([5]) coming from a sum over the colors of the c, and c, making up the J/^- To find the color 
factors in the gluon induced process it is necessary to evaluate the lowest order process shown in 
Fig. [2] The color factors for this process are, as already found in [2], 

^^T. IV («'<«) . (^) " ^ (12) 



where we have used the large- A'^c hmit in the right hand side of (12). The (^Cp/{N^ — 1)) factors 
go into making up part of the two factors of that come from the graphs. Explicit calculation 



confirms that the remaining factor, after taking out the factor in (^9| and the factor linear in Q 



when expanding (10) is just the factor l/2Nc on the right hand side of (12). 



Putting all this together gives 



'gA^J/^liX 



10 J "t" JO 

X / di Z^^" exp 

^0 



(fr' 



47r In I Att 



4r(5) 



1 

16 



T{b) 8 



T{b) 



(13) 



with 



<I>(r, z) 



:eN, 



(14) 



where, finally, in (14) we have introduced the replacement |e — ?■ g. Doing the integral over ^ and 
using ([1]) we get 



4r • r' 
(r + r')2 













47r 







-%{r2+r'2) 



(15) 
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III. J/V- CROSS SECTION IN AA COLLISIONS 

Generalization of the result of the previous section to nucleus-nucleus collisions is achieved by 
letting the initial gluon be emitted from either nucleus and taking into account cc dipole scattering 
in both nuclei. The scattering amplitudes and the saturation scales for the two nuclei depend on 
their respective impact parameters hi and b2- To make our notations more compact we will not 
indicate the impact parameter dependence explicitly. Introducing the relative impact parameter 
B = b]^ — b2 and using the relation 

^'^\iG{xi,a'^)=[d%^ (16) 



ACf 

we can write the cross section as 



da^,^,^j,,x .J^^Jdhj^' ^, I ^^^^ 2T^,^,^Mr_, r_') , (17) 



dy(Ph(PB 

where 

-,2 ri2 



Ta a ^lYii r') - 4r-r^ / _^(Q2^+Q2^)(r_r')2 _ -i{Q2^+Q2^)(r2+r'2) 



(18) 



Expanding (18) at small Q^^i we recover Eq. (15) 



The first few terms in the expansion of (18) in nuclear density read 

Ta,A,^jx{x, r') ^ ^^QliQl^ 4r • E - ^(Q^i + Q?2)(3r2 ^ 3^/2 _ 3^ . j (19) 

Averaging over the relative angle between r and r' yields 

(rAM2^Jx(r,r')> « (QfiQ!2 + Q'iQ'2) (20) 

yJig TV J. O 

This is the leading contribution to the J/tp production; it is easily seen that it breaks the factor- 



ization. We believe that (18) is a reasonable starting point for phenomenology of J/ip production 



in AA collisions. Nevertheless a better theoretical understanding of the AA production amplitude 
r^,^2^j/^x is desirable. 



IV. RAPIDITY AND ENERGY DEPENDENCE 



Eqs. (17), (18) can be readily generalized to include quantum evolution effects. To that end we 
recall that the initial condition for the BK [5l [10] evolution equation is given by the Glauber-Mueller 
formula for the forward dipole-nucleus quark dipole elastic scattering amplitude yO] 



NF{r,h,yo) = l-e-"s^-'^"^yo) , (21) 
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where subscript F indicates the fundamental representation. Evolution of the gluon dipole scat- 
tering amplitude (adjoint representation) obeys the equation 



Na{t, b, y) = 2Nf{t, b, y) - N^l, b, y) 



and its initial condition is 



NAiL,h,yo) 



1 _e-3E'Q?feo) 



Accordingly, we can incorporate evolution effects in ( 18 ) by the following replacements 



e 16^ n 



l-NA{l/2,h,y) 



(22) 



(23) 



(24) 
(25) 



Omitting the impact parameter dependence as before, we thus obtain 
M _ 8N, Q\Ql^ 4r • r' 



rAiA2^Jx(E,E ) 



a.TT^ Q\^Ql^ (r + rO 



^ { [l - Nf{{j_ - i)l2, y)\ [l - ivf ((r - r')/2, -y) 
\-f^f(j_,y)\ [l-iv(|)(r',y)] [l - Ar(2)(r, -y)] _ Ar(2)(,', | 



(26) 



V. NUMERICAL CALCULATIONS 



The experimental data is expressed in terms of the nuclear modification factor (NMF) defined 



as 



(27) 



dy 



where S stands for the overall area of two colliding nuclei. Since the mechanism of J/ip production in 
pp collisions remains elusive, we follow our approach in the previous publications and approximate 



da. 



pp^J/tpX 



C 



da 



AA-^J/ipX 



(28) 



A=l 



dy dy 

with C = const. We fix the constant to provide the best description of the pp and dA data. It is 
reassuring that the numerical calculations described in the next section indicate that C is close to 
unity. 

The results of our calculations are exhibited in Fig. [3] and Fig. [4j we have used two different 
models for the dipole scattering amplitude: DHJ [T3] and bCGC models (see Appendix |A] 
for the description of these models). Comparison of the results of the two models gives an idea 
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about the model dependence of the numerical results. We observe a reasonable agreement with 
the experimental data on J /ijj production in dA collisions. 

Concerning the J/ij: production in AA collisions all models underestimate the suppression at 
RHIC both at mid-rapidity and in the forward rapidity. Moreover, it appears that the gluon 
saturation effects on NMF show very little rapidity dependence at RHIC which contradicts the 
experimental data. We also find that there is almost no change between the NMF at LHC = 
2.76 TeV and 5.5 TeV. We note that our calculation overestimates the NMF at = 2.76 TeV. 

1.2 : 
1.0- 
o.i 



10 



12 



14 



16 



N„. 




(a) (b) 

FIG. 3: Nuclear modification factor vs iVpart in (a) dAu and (b) AA collisions using the DHJ model [T5] . 
Band 'a' (green) represents rapidity y = —1.7 at -y/s = 200 GeV, 'b' (blue): y — — 200 GeV, 'c' 

(red): y = 1.7, y^s = 200 GeV, 'd' (brown): y = 3.25, = 2.76 TeV, 'e' (cyan): y = 0, = 5.5 TeV. 
m = 1.5 GeV, C = 1. Experimental data [TBHT^ is represented by (blue) circles in 'b', by (red) squares in 
'c' and by (brown) triangles in 'd'. (Color online). 
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VI. DISCUSSION AND CONCLUSIONS 



Our calculations indicate that the nuclear modification of J/tp production in dA collisions at 
RHIC is dominated by the cold nuclear matter effects. It would be important to study J/ip 
production in pA collisions at LHC; Fig. [3] and Fig. [4] provide our predictions. In contrast, the cold 
nuclear matter effects alone cannot provide neither quantitative nor even a qualitative description 
of the AA data. Additional mechanisms beyond the initial state effects are required to explain 
the experimental data. It is remarkable that at RHIC these additional mechanisms must provide 
extra suppression of the NMF, perhaps via the Matsui-Satz color screening mechanism or the 
gluon-induced dissociation [24^ I25j. whereas at LHC they must produce enhancement. 

Our successful description of the J/ ij: NMF in pA collisions with the normalization factor C = 1 



in ( 28 ) may be an evidence that the J/tp production mechanism in pp collisions is similar to that 
in pA implying that it is perhaps dominated by the higher twist effects. 

To summarize, we derived the formulae for the cross sections of J/ijj production in pA and AA 
collisions taking into account the gluon saturation/color glass condensate effects. Our numerical 
results provide an estimate of the color nuclear matter effects on J/ip production in heavy-ion 
collisions. 
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Appendix A: Models of the dipole scattering amplitude 

We performed numerical calculations using two models of the dipole scattering amplitude: DHJ 
|13j and bCGC [9j models. The DHJ model is an improvement of the KKT model[T31 EI] that 
takes into account the change in the anomalous dimension of the gluon distribution function due 
to the presence of the saturation boundary [15] and takes into account some higher order effect. 
It successfully describes the single inclusive hadron production in dA collisions in the relevant 
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kinematic region. In this model, the dipole scattering amphtude is parameterized as follows 

1 ' 2^2^7 



iVA(r,0,y) = l-exp|-- (r^Q^)'j . (Al) 
The gluon saturation scale is given by 

= A^/^ = 0.13 GeV^ iVcoU • (A2) 

where the parameters A = 0.6 GeV and A = 0.3 are fixed by DIS data |22| . 

ln(mVQ?) ^ 

7 = 7, + 1-7, ^ '-^ — A3 



where 7^ = 0.628 is implied by theoretical arguments [15] and d = 1.2 is fixed by fitting to the 
hadron production data in dA collisions at RHIC. Y = ln(l/x), with x = me~^/-y/s. The quark 
dipole scattering amplitude is given by 



Nf{t. 0, y) = 1 - Vl-A^A(r,0,y) (A4) 



which follows from ( 22 ) . 



We used the bCGC model [9j with a modification: we treat the nuclei and proton profiles as 
step-functions; the saturation scales are assumed to scale with A as Ql on A^^^. The advantage of 
this model - besides its compliance with the known analytical approximations to the BK equation 
|23j - is that its parameters are fitted to the low x DIS data. The explicit form of the scattering 
amplitude is given by 

iVF(r,0,y)= <^ V 4 ; ' (A5) 

[ l-exp[-aln2(6rQ,)], rQ, > 2 , 

where is the the quark saturation scale related to the gluon saturation scale - which we 
have called simply the 'saturation scale' throughout the paper - by Q2 = (4/9) Q2. Its functional 
form is 

where s is the square of the center-of-mass energy and y is rapidity with respect to the central 
rapidity. The anomalous dimension is 

c\{\n.^Js + y) \rQsJ 



11 



Parameters 7^ = 0.628 and c = 9.9 follow from the BFKL dynamics [23], while A/q = 0.7 and 
A = 0.28 are fitted to the DIS data. Constants a and b are uniquely fixed from by the requirement 
of continuity of the amplitude and its first derivative. 
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